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Abstract. The evolution of masses and sizes of passive (early-type) galaxies with redshift pro- 
vides ideal constraints to galaxy formation models. These parameters can in principle be ob- 
tained for large galaxy samples from multi-band photometry alone. However the accuracy of 
photometric masses is limited by the non-universality of the IMF. Galaxy sizes can be biased 
at high redshift due to the inferior quality of the imaging data. Both problems can be avoided 
using galaxy dynamics, and in particular by measuring the galaxies stellar velocity dispersion. 
Here we provide an overview of the efforts in this direction. 



Keywords, galaxies: elliptical and lenticular, cD - galaxies: evolution 
galaxies: kinematics and dynamics - galaxies: structure. 



galaxies: formation 



1. Introduction 

In the era of precision cos mology we think we c an accurately predict the distribution of 
dark matter in the Universe (jSpringel et al ] l2005l) . Dark matter is thought to be the main 
driver for the hierarchical assembly of galaxy stellar masses. However the mechanism by 
which the gas sinks into the dark matter potential wells and form stars is still largely 
unknown, due to the complex and difficult to model effect of baryonic physics. For this 
reason our understanding of galaxy formation must be driven by observations. 
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2. Uncertainties of galaxy stellar masses 

A first key observable to compare with galaxy formation models is the galaxy stellar 
mass, which is expected to monotonically growth with time during galaxy mergers. In 
recent years confidence started to growth about our ability to measure stellar masse s 
with a few tens of percent accuracy (e.g. iGallazzi fc Bell 120091 : Mar aston et a l. 2010), 
up to redshift z ~ 2 and beyond, thanks to the availa bility of multi-bands p hotometric 



surveys with the Hubble Space Telescope (H TS) (e.g. iGiavalisco et alll2004l) combined 



with detailed stellar population models (e.g . Bruzual fc Chariot! 120031 : iMarastonl 120051 : 



with detailed stellar population models (e.g . rJruzi 
Vazdekis et al.ll20ld IConroy fc van DokkumlkoA 



Mass determinations using stellar population models depend on the assumption of 
the stellar initial mass function (IMF). Unt il recently it ap pe ared sensible to assume the 
IMF is universal in different galaxies, and of lKroupal (|200lh or IChabrierl (|2003h type. This 
assumption was motivat ed by the observed IMF universa lity in different environments 



within our own Galaxy (jKroupa 2002; B astian et al.l l2010) 



Recent results from spectral absorption featur es (e.g. van Dokkum fc Conroy 2010h 
and gravitational lensing (e.g. lAuger et al. 2010l) found indications for a hea vier IMF 
in ma ssive early- type galaxies (ETGs), more consistent with an IMF with the ISalpeter 



(|l955l) ogarithmic slope or even heavier, placing the IMF universality into question. The 
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extensive dynamical modelling effort of ICappellari et al.1 (|2012al) discovered a systematic 
trend in the IMF of ETGs. They found a mass normalization for the IMF varying between 
values consistent with a Kroupa/Chabrier IMF to heavier than a Salpeter IMF (assuming 
the population models are correct) as a function of either the mass-to-light ratio (M/L) 
or the galaxy stellar velocity dispersion a. This trend implies potential systematic errors 
of up to a factor 2-3 in the galaxies stellar masses derived using stellar population models, 
even from optimal data. 



3. Uncertainties of galaxy sizes 

A second observable that is often used to constrain galaxy formation scenarios, in 
combination with galaxy mass, is the galaxy size, commonly parametrized by the pro- 
jected half-light radius R c . The evolution of galaxy size with its mass depends sensitively 
on the galaxy assembly mechanism (see Ciotti 20091 for a comprehensive review). R c 
increases nearly proportio nally to mass in the ca se of a galaxy assembly process via gas- 
poor major (1:1) mergers (jHernquist et al. 1993), due to energy con servation, while R e 
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increases more rapidly dur ing gas-poor minor (1:3 or less) merging (jNipoti et al. 
Bovlan-Kolchin et al. 2006( ). R e decreases during gas-rich mergers o r cold accretion, when 
gas accumulates towards the centre (e.g. iMihos &; Hernquisdll994 ) . 

Massive (stellar mass M* > 10 11 M©) passive galaxies at z ~ 2 were found to have 
much smaller R e than their local counterparts of the same mass (LDaddi et al. 2005; 



di Serego Alighieri et alj2005l:lTruiillo et al.ll2006l;lLonghetti e t alJl2007t lToft et al.l l2007; 



Cimatti et all 120081; Ivan Dokkum et al.1 120081) ■ This suggested a scenario in which the 



prog enitors of today's massive ETGs assembl e most of their mass v ia minor dry mergers 
(e.g. Naab et al. 2009t Bezanson et al. 2009t Hopkins et al. 20091 ) . to account for the 
rapid size increase. 

Unfortunately the half-light radius R e is an observationally ill-defined quantity. By 
definition measuring R c requires an extrapolation of the galaxy surface brightness profile 
to infinite radii, to estimate the galaxy total stellar light. This seemed a sensible approach 
a few decades ago, when the technique was used to quantify the sizes of elliptical galax- 
ies, w hich were thought to be all well described by homologous i? 1//4 Ide Vaucouleurs 
( 19481) su r face-b rightness profiles. But we now know that ETGs are best represented 
by ISersic (Il968l) profiles, with concentration increasing with galaxy mass (|Caon et al.1 
1993]). Moreover the ATLAS 3D volume-limited (M* > 6 x 10 9 M ) survey of nearby 



ETGs (|Cappellari et al.l 1201 lal) has revealed that the majority of the s ystems are more 
closely related to spiral galaxies than to genuine spheroidal ellipticals ( Cappellari et al.1 
2011bl) . In fact as much as 2 /3 of the galaxies classified as ellipticals are dominated by 



stellar disks, which are often mi ssed by the photometry, but are clearly revea led by their 
integral-field stellar kinematics (jKrainovic et al. I bOllHEmsellem et aijlioil . 

The complex and multi-component nature of ETGs makes it impossible to know how 
the outer unobservable surface brightness profile should be extrapolated. Different au- 
thors use photometry of different quality and make different choices for the profile ex- 
trapolation, as well as for its extraction. For these reason, even from very good quality 
photometry of nearby galaxies, R e is often not known with an accuracy better than 
about a factor 2x, with error s dominated by s ystematics effects. This is illustrated in 
the left panel of Fig. Q] (from lChen et al.ll2010l ). which provides a comparison between 
independent determinations of R e for well-studied ETGs in the Virgo cluster. 

Uncertainties and systematic biases in R e are likely to be more severe when galaxies 
at high redshift are compared to local ones. Differences in the measurement approach 
combine with (i) differences in the restframe wavelength of the observations, (ii) variations 
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Figure 1. Left: Uncertainty of R c for galaxie s in the V i rgo c luster. Top Panel: Different 
symbols compare the SDSS half-light radii from Ch en et al.l l|2010f ) with independent determi- 
nations of the same quantity from the literature. Bottom Panel: Gray dots show the fractional 
error in the effective radius measurements, plotted as a function of the SDSS effective radii; the 
mean and standard deviation for individual galaxies are plotted as the black dots with error 
bars. Differences in Re of a factor 2x are common, especially for the largest galaxies, (taken 
from [Chen et al1l2010f h Right: Comparing images of nearby and distant ETGs. Top 
Panels: true colour images of n earby ETGs from the ATLAS 3D sample as imaged by the SDSS 
survey l|Abazaiian et all [2009). Bottom Panels: true colour images of galaxies at z = 1 — 2.5 
from the CANDLES survey (taken from ISzomoru et al.ll2012T) 



in the depth of the photometry (e.g. due to cosmological surface-brightness dimming) and 
(iii) inferior spatial sampling of the imaging of distant galaxies (see the right panel of 
Fig. [JJ, (iv) possible colour gradients and nuclear AGN activity, which are generally 
stronger in the early Universe. 



4. Local benchmark for dynamical scaling relations 

An alternative to measuring galaxy sizes to study the evolution of galaxy densities 
consists of extracting the galaxies velocity dispersion a. This is more difficult to obtain 
than sizes, as it requires spectroscopic data (resolution R w 1000) rather than imaging 
alone. However a does not suffer from the biases of R e and is closely related to it via the 
scalar virial equation (e.g. Cappellari et alJfefJOd) 



yn 



5.0 R e 



(4.1) 



Here Md yn = 2 x M 1 / 2 ~ M*, where M 1 / 2 is the mass within a sphere enclosing half of the 
total galaxy light , and the last approxima tion is due to the fact that Mi/ 2 is dominated by 
the stellar mass (ICappellari et al.ll2012ch. The factor 5.0 was calibrated for R e measured 
in the classi c way dBurstein et al.ll987t de Vaucouleurs et al.lll991 ; Jorgensen et aljl995l) 
using fixed Ide Vaucouleurs! (|l948l) R ^ growth curves to extrapolate the outer profiles 
and using typical photometry of nearby galaxies. 
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Figure 2. The Mass Plane and it projections. The top two panels show the (MjAM,o" e ) 
and (MjAM,i?r ax ) coordinates. Overlaid are lines of constant cr e = 50,100,200,300,400,500 
km s _1 (dashed blue), constant _R™ ax = 0.1, 1, 10, 100 kpc (dot-dashed red) and constant 
E e = 10 s , 10 9 , 10 10 , 10 11 M Q kpc" 2 (dotted black) predicted by the virial equation (|4~T)) . The 
observed (Mjam, o" e , _R" lax ) points follow the relation so closely that the coordinates provide a 
unique mapping on these diagram and one can reliably infer characteristics of the galaxies from 
any individual projection. In each panel the galaxies are coloured according to the (LOESS 
smoothed) log(M/L)jAM values, as shown in the colour bar. Moreover in all panels the thick 
red line shows the same ZOE relation projected via equation (14.1 [I. The g reen line is the M — R c 
relation for late spiral galaxies (equation 3 from ICappellari et alJl2011ah. which ap proximately 
defines the boundary where ETGs disappear. (Taken from ICappellari et al.ll2012bf ) 



To study the evolution of galaxy parameters with redshift, a reliable local benchmark 
is essential. In Fig. [2] we show the (M, a) and (M,R e ) (and other projections) of the 
mass plane (M, R e , a) of ETGs obtained by the AT LAS 3D project via detailed dynamical 
modelling of 260 galaxies (jCappellari et al. I l2012bl) . The study found that, when accurate 
and unbiased masses are used, ETGs satisfy equation (|4. 1[) quite accurately. Different two- 
dimensional projections of the plane contain the same amount of information except for 
a coordinate transformation. This confirms in particular that one can reliably estimate 
a from the knowledge of M and R c using equation (|4.1j) . The study also found that the 
zone of exclusion (ZOE) defined by local ETGs in the (M, R e ) plane shows a clear break 
at a characteristic mass Md yn ~ 3 x 10 10 M©, with a corresponding break in the ZOE in 
the (M, a) plane and other projections. 
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Figure 3. Mass, size and a distributions at ; 

terminations (black fi l led circles) for ETGs with 1.4 < z < 2.5 (from [Cenarro fc Truiillc 



2. Top Left: The observed (M,q) de 



Cappellari et al.ll2009l; Ivan Dokkum et ai1l2009l ; Ivan de Sande et al.ll2011l ; iNewman et al 



2009 



2010 



Onodera et al.ll2012l ; lToft et al.ll2012l ) is overlaid to the locus of local ETGs from ATLAS 3D (grey 



region). The red and green thick lines are the same as in Fig. [2] Top Right: as i n the previous 
pane l , for the distribution of (M, R e ) determinations (blue filled circles) (from ICassata et al.l 
2011; Szomoru ct al. 2012). Bottom Left: The (M, a) distribution predicted using equation (|4.1|l 
from the values in the top-right panel. Bottom Right: as in the previous panel, but using 2R C to 
estimate a instead of the observed R e . 



Indications for a change of slope in the Faber fc Jackson! 



. .. , - 19761) relation between galaxy 

luminosity and a ( Davies et al. 1983 ) or in the iKormendv ( 1977 ) relation between lumi- 
nosity and surface brightness ( Binggeli et al. 1984) were found some time ago. However 
this is the first time the break is presented using dynamical masses instead of luminosi- 
ties, and it is shown to accurately represent two projections of the same ZOE on the 
galaxies mass plane. This break is interpreted as evidence for a change in the galaxy ac- 
cretion mechanism between bulge growth (in-situ star formation) and gas-poor merging 
(external accretion). 



5. Importance of o determinations at high-z 

In the top-left panel of Fig. [3] we sh ow the available determinations of a for galaxies 
in the redshift interval 14 < z < 2.5 dCenarro fc Truiillol 120091: ICappellari et all 120091: 
van Dokkum et "all 120091: Ivan de Sande et al.ll2011t INewman et al.ll2012t I Onodera et al 



20121: iToft et al.l 120121) and compare it with the location defined by nearby ATLAS 3D 
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ETGs. There is an indication for the high-redshift a determinations to lie generally 
above the locus of local galaxies with the same mass, as found earlier, and as expected 
for more compact systems. However, with the exception of one galaxy, the differences are 
relatively modest and are mostly consistent with local galaxies within the measurement 
errors. 

In the top-righ t panel of Fig, p we show state-of-the-art size determinations of R e for 
passive galaxies ([Cassata et al. zomoru et al.ll2012l ). in the same redshift range. 

They were obtained from deep HST/WFC3 observations in the GOODS and CANDLES 
fields. The photometric samples are much larger than the one with a determinations, 
and in this case the difference between the locus of local ETGs and the z ~ 2 ones is 
more dramatic, with the high-z galaxies having much smaller sizes than local ETGs. 

In the bottom-left panel of Fig. [3]the observed (M, R c ) distribution of the z ~ 2 ETGs 
is converted into the (M, a) plane via equation (|4.1|) . This shows that the observed dis- 
tribution of galaxy sizes implies more extreme a values than the ones that have been 
found so far. The expected a may even be larger than predicted here, due to the fact 
that high-z galaxies not only are observed to be smaller, but a lso tend to have steeper 



profiles than their local counterparts (e.g. ISzomoru et al.l 120121) . One way to bring the 



high-z predicted a values into better qualitative agreement with the few available deter- 
minations, would be to assume the R e values at z ~ 2 are underestimated by about a 
factor of 2x (bottom-right panel of Fig. |3|). 

However the comparison between the currently available a and R c determinations at 
z ~ 2 is not sufficient to conclude that the high-z determinations of R c are underesti- 
mated. The sample of galaxies with available a is still too small to draw robust con- 
clusions, and it is selected in a very different way than the more complete photometric 
samples. Moreover the galaxies with available a generally satisfy equation (|4.1|) reason- 
ably well. What the comparison does emphasizes is the importance of obtaining more 
good-quality a values of high-z galaxies (see e.g. van de Sande in these proceedings), 
until a fully consistent picture emerges. 

The possibility of biases in the R e determinations of high-z galaxies also emphasizes 
the danger of using the local virial equation (|4.1j) for quantitative estimates of M/L, e.g. 
to study the evolution of the Fundamental Plane or to explore possible variations of the 
IMF with redshift. There is a clean solution to the problem. It consists of using detailed 
axisymmetric dynamical models of high-z galaxies (lyan der Marel fc: van Dokkuml 12007 ; 
van der Wei fc van der Marellliooj ICappellari et al.ll2009f )7 which describe in detail the 



shape and surface brightness profiles of the galaxies. Contrary to the dynamical M/L 
inferred from the virial relation, the M/L derived via dynamical models is virtually 
insensitive to even extreme changes in the surface brightness distribution outside the 
region where stellar kinematics is available. Alth ough the correctness of this fact is best 
demonstrated with numerical experiments (e.g. Cappellari et al. 20091) . the reason for 



this robustness is similar to the one for which, in the spherical limit, the knowledge of 
the orbit of a single star is sufficient to infer the mass enclosed inside its orbit, irrespective 
of the large-scale mass distribution. 

A demonstration of this modelling approach is given in the left panel of Fig. 21 
which shows a first att empt at measuring the IMF normalization for galaxies at z ~ 2 
( Cappellari et a02009() . The study used dynamical models to reproduce a determinations 
and concluded that some of the high-z ETGs are required to have a mass normalization 
of the IMF as low as Kroupa/Chabrier, like local ETGs. However the study allows for 
the possibility of some of the studies galaxies to have a heavier Salpeter-like mass nor- 
malization. This is not inconsistent with local findings, which indicate the IMF should 
vary between Kroupa to Salpeter in the interval a w 80 — 260 km s _1 (right panel of 
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Figure 4. Constraining the IMF at z ~ 2. Le/t Panel: dynamical versus population masses 
determined using a Chabrier IMF. Blue diamonds are galaxies with individually measured a, 
while red filled circles are galaxies for which a from a stacked spectrum was assumed. The green 
solid line is the one-to-one relation. The dashed line is the one-to-one relation if the IMF was 
of Salpeter type. The latter is inconsistent with some of the determinations, i ndicating that 
some of the object require a Chabrier-like IMF (taken from Cappellari et al.ll2009f) . Right Panel: 
Location of the galaxies in the left panel in the (M, a) plane. The grey area and thick red 
and green lines are as in Fig. [2] The l abels show the IMF nor malization expected at different 
locations on this plane for local ETGs l|Cappellari "et~alll2012bl ). 



Fig. H iCappell ari et al. 2012b). But the current data and samples are still not sufficient 
to draw strong conclusions on the IMF variation at high-z. 

The near future looks bright for dynamical studies of galaxies at z ~ 2, thanks to 
the availability of multi-object Spectrometers like MOIRCS on the Subaru telescope, the 
upcoming MOSFIRE on the Keck telescope and KMOS on the Very Large Telescope. 
These are optimized for multiplexed near-infrared observations of distant galaxies and 
should provide in a few years numerous a determinations at high rcdshift ETGs. 
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